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ABSTRACT 

We investigate the alignment processes of spinning black holes and their surrounding warped accretion disks 
in a frame of two different types of feeding at the outer boundaries. We consider (1) fixed flows in which gas is 
continually fed with a preferred angular momentum, and (2) free flows in which there is no gas supply and the 
disks diffuse freely at their outer edges. As expected, we find that for the cases of fixed flows the black hole 
disk systems always end up aligning on timescales of several 10 6 yr, irrespective of the initial inclinations. If 
the initial inclination angles are larger than ir/2, the black hole accretion transits from retrograde to prograde 
fashion, and the accreted mass onto the black holes during these two phases is comparable. On the other hand, 
for the cases of free flows, both alignments and anti- alignments can occur, depending on the initial inclinations 
and the ratios of the angular momentum of the disks to that of the black holes. In such cases, the disks will 
be consumed within timescales of 10 6 yr by black holes accreting at the Eddington limit. We propose that 
there is a close connection between the black hole spin and the lifetime for which the feeding persists, which 
determines the observable episodic lifetimes of active galactic nuclei. We conclude that careful inclusion of the 
disk feeding at the outer boundaries is crucial for modeling the evolution of the black hole spin. 
Subject headings: accretion, accretion disks — black hole physics — galaxies: active 



1. INTRODUCTION 

The current hierarchical framework of galaxy formation 
and evolution predicts that, in addition to secular evolution, 
repeat galaxy mergers trigger activities of supermassive black 
holes (SMBH s) located at the galactic center and hence their 
growth (e.g., iBenson & Bowed |2010|) . As such, fuel feed- 
ings channeled onto SM BHs most likely proceed at episodic 
and random phases (e.g. , lMa rtini 2004; King & Pringle 20061: 
IWang et alJ [2006L [20081 [20091: This picture 

is further reinforced by various lines of observations, such 



their host galaxies ( 


Kinnev et al. 2000; Gallimore et al. 2006; 


iMunoz Marin et al. 


2007; Shenetal. 201(3), and as well as 



erences therein). Previous studies on the cosmological evolu- 
tion^ the radiative efficiency (77) of accretion according to the 
ISoltanl 's (1982) argument and formulated by the /^-equation, 
found that 77 decreases with cosmic time since redshift z ~ 2, 
strongly implying that random accretion takes pl ace upon the 
SMBHs ( Wa ng et al.ll2009l: iLi et aUl20TIL 12012 : see also the 
discussion of IZhang et al.ll2012l) . A similar inference has re- 
cently been highlighted from the semi-a n alytic al modeling on 
SMBH demography bv IVolonteri et al.1 (120121) . In this con- 
text, warped accretion disks are ubiquitous in AGNs. 

A misaligned accretion disk around the central spinning 
black hole undergoes Lense-Thirring precession arising from 
the gravitomagnetic effects, of w hich the rate falls off rapidly 
with the radius to the hole (e.g., Hartle 2 0031) . The presence 
of viscosity combined with such differential precession will 
induce the inner portion of the disk to align or anti-align its or- 
bital angular momentum with the spin axis of the hole out to a 
transition radius, beyond which the disk retains its initial incli- 
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nation (Bardeen-Petterson effect: iBardeen & Pettersonl 19751) . 
In the meantime, associated with disk precession, the black 
hole suffers an equal and opposite gravitomagnetic torque 
that causes it to precess as well. Under mutual interac- 
tion and warp propagation in the disk, this composite sys- 
tem tends to restore full axisymmetry with the entire disk 
ultim ately aligned or anti-aligned with the hole dKing et al.l 
120051) . The characteristic timescale and the final state of the 
system (alignment or anti- alignment) depend on the black 
hole mass and spin and the properties of the accret ion disk. 

Bas ed on the angular momentum conservation, iKing et ail 
(120051) demonstrated that the disk will end up anti-aligned 
with the hole if the angular momentum of the hole (/h) dom- 
inates over that of the disk (J&). This refreshes previous an- 
alytical studies on the grounds that the disk has an infinite 
extension and is continuously fed (IScheuer & Feilerl fl996t 
iMartin et"aLll200 7). and has been incorporated into the black 
hole spin modeling in the context of hierarchical cosmology 
dKing et al.1 [20081 : lLagos eTall [20091 : iFanidakis et all [20111: 
iBarau sse 2012). However, while Jh is acc urately defi n ed, th e 
meaning of J& is not as well defined in IKing et ail (120051) . 
leading to its interpretation varying among subsequent stud- 
ies. 

It is commonplace that, in the outer region, the accre- 
tion disk likely becomes self-gravitating, fragments into 
clumps, an d probably manifests itself as a star-formi ng re- 
gion (e.g., Ea czvnski 197 8]; [Shj osman & Begelman 1987; 
iCollin & Zahnl 119991: iGoodmanl 120031: IWang et al.1 12010b . 
The feedback from the stellar radiation or the resultant su- 
pernova explosions embedded in the disk acts to excite 
turbulenc e that continues to deprive gas of angular mo- 
mentum dWada & Normanl 120021: iKawakatu & Wadal 120081: 
iKumar & John son 20 lot I Wang et al.ll2010l) . In this case, it is 
expected that gas will be continually fueled into the accretion 
disk from the outer region (e.g., circumnuclear disk or torus) 
for the duration of the AGN's lifetime. The practical impli- 
cations are twofold: an outer edge of the disk limited by its 
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own self-gravity exists so that the angular momentum of the 
disk can be properly defined, and there is a gas supply with 
specified angular momentum at this outer edge. 

With the above inference, a comprehensive study of the 
alignments/anti-alignments between black holes and their 
warped accretion disks is essential, in particular considering 
the increased focus on modeling the black hole spin across 
cosmic time. In this work, we revisit this issue by numer- 
ically solving the evolution equation set of black holes and 
accretion disks in a practical sense. The paper is organized as 
follows: Section 2 describes the evolution equation set for a 
black hole and its warped accretion disk. Section 3 surveys 
the properties of warped accretion disks from the theoretical 
point of view. Section 4 shows the numerical scheme and sets 
up the boundary conditions and initial conditions for solving 
the evolution equation set. The results are then presented in 
Section 5. The implications of our results on SMBH spin evo- 
lution and the uncertainties in our calculations are discussed 
in Section 6. Conclusions are summarized in Section 7. 

2. BASIC EVOLUTION EQUATIONS 

The dynamics of warped disks hav e been studied exten- 
sively through theoretical analysis (e.g.,|Pa paloizou & Prmgl 



through theoretical analysis (e.g., rapaloizou & Fnnsle 
lPringld[T99l iPaoaloizou & Linl 119951: lOgilyid I1999L 
lLubow et al.l |2002|) and numerical simulatio n s (e.g . 
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iNelson & Papaloiz oul [T999i 



2000; 



2012; 



Dotti et alJl2010b lLodato & Price 2010; Nixo n & King _ 
Nixon et al.ll2012l:lMcKinnev et alJl2013h . Adopting the stan 
dard a-prescription of disk viscosity, warp propagation in the 
disk is characterized into two regimes depending on the rela- 
tive importance of pressure forces and viscous forces. For a 
nearly inviscid or sufficiently thick disk such that the dimen- 
sionless viscosity coefficient a is smaller than the disk aspect 
ratio h = H/R, the warp propagates as bending waves at a ve - 
locity related to the speed of sound ( Papaloizou & Lin|[l995|) : 
whereas in the regime a > h, the warp evolves in a diffusive 
fashion. Here H is the semithickness of the disk. 

We focus our attention on a thin and viscous accretion disk 
with a Keplerian rotati on (a > h), the most likely case for ac- 
cretion disks in AGNs (Pringle 1999). The evolution of such 
warped disks can be completely descri bed by the ang ular mo- 
mentum density of the disk L(R,t) as (10 gilvidl 1 999|) 
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where L = Y*R 2 Q, £ = h/L denotes the direction vector, E is 
the mass surface density, ft is the angular velocity, Tlt is 
the gravitomagnetic torque responsible for the Lense-Thirring 
precession, and Vr is the radial velocity given by 



3 d(y x L) 

Vr = uiR 

L dR 



dl 



dR 



(2) 



Here the internal viscous torque of the disk fluid has been 
distinguished into three kinds associated with three effective 
visco sities, respectively dBardeen & Pettersonlll975l: lOgilvid 
119991) . The first viscosity, v\, governs the usual azimuthal 
shear due to the differential rotation; the second viscosity, 
V2, governs the vertical shear that acts to diffuse the warp 
throughout the disk; and the third viscosity, ^3, governs a 
torque component that tends to create the disk ring process 



if it is misaligned with its neighbors. The later two viscosities 
disappear in a flat disk. It is worth pointing out that Equa tion 
(Q} remains valid for large amplitude warp (Ogilvie 1999)). 

By analyzing the inte r nal no n-linear dynamics of warped 
accretion disks, lOgilvid (Il999|) formulated the effective vis- 
cosities and showed that they generally depend on the ampli- 
tude of the warp. For the sake of simplicity, we make use of 
the first-order approximations as follows 



v x = ah 2 R 2 n, 



(3) 



V2 — v\j 2a 2 and ^3 = 3 v\ /8a. As a result, the warp propagates 
in the disk on a timescale of an order of 



R 2 r> 7 

t w = — ~ 2a — =2a t v , 
where t v is the usual viscous timescale. 



(4) 



In the weak field limit, the gravitomagnetic torque du e to 



frame dragging has a form of dBardeen & Pettersonl 



q ue du e 
n ll975h 



Tlt = 



2G J h xL 

c 2 R 3 ' 



(5) 



where Jh is the angular momentum of the hole with its mod- 
ulus /h = aGM 2 /c, a is the spin parameter, M m is the mass of 
the black hole, G is the gravity constant, and c is the speed 
of light. Due to the mass accretion and the gravitomagnetic 
torque exerted by the warped accretion disk, the angular mo- 
mentum of the black hole evolves with time as 



dJh 

dt 



:M in j in -2^ / T LT RdR, 



(6) 



where M- m is the mass accretion rate through the inner edge of 
the disk R m , which is determined by the marginal stable orbit, 
and ji n is the specific angular momentum carried by matter at 
the inner edge. 

To study the roles of disk feeding, we consider two classes 
of accretion models to mimic the realistic accretion disks, as 
schematically shown in Figure [T] In the left panel, an outer 
clumpy region composite of clouds and/or stars continues to 
channel gas to the accretion disk from the outer environment 
(e.g., circumnuclear disk or torus). As stated above, the feed- 
back from stellar radiation and supernova explosions are the 
potential mechanism responsible for exciting turbulence and 
viscosity, a cting to transport the ang ular momentum of accret- 
ing matter (IWada & Normanll2002t IKumar & Johnson! l2010t 
Wang et al. 2010|). We denote this case as "fixed flows" be- 
cause the angular momentum of gas supply at the outer edge 
has a preferred direction. By contrast, the right panel of Fig- 
ure [T] shows the other class of outer boundary conditions, in 
which there is no gas supply and the disk diffuses freely at its 
outer edge. We therefore denote this case as "free flows". Be- 
fore numerically solving the above equations, we endeavor to 
explore the properties of warped accretion disks from a theo- 
retical point of view in the next section. 

3. WARPED ACCRETION DISKS 

We start by specifying t he accretion disk model used in thi s 
paper generally following ICollin-Souffrin & Dumontl (119901) . 
To simplify the theoretical analysis and numerical calcula- 
tions below, we parameterize accretion disks with a Keple- 
rian rotation of ft = y/GM 9 /R 3 , and a constant aspect ra- 
tio of h = H /R. The mass surface density of the disk obeys 
E = So^ _/? , where So is a coefficient determined by set- 
ting the total disk mass at a fraction of the black hole mass 
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Fig. 1. — Schematic illustration of two classes of accretion models. Left: fixed flows in which gas is continually fed with a preferred angular momentum 
through the outer region composite of cloud and stars. Right: freeflows in which there is no gas supply and the disk diffuses freely at its outer edge, leading to 
an outflow of angular momentum. The central black hole spin is oriented with an inclination to the accretion disk. The outer edge of the accretion disk is limited 
by self-gravity of the disk and the inner edge is determined by the marginal stable orbit. 



Md = f m M m , and p = 1/2 is an index to enforce the mass ac- 
cretion rate of M « 37r^iE to a constant at a large radius as 
for a flat disk. The outer edge of the disk R out is located where 
self-gravitating becomes important, whose value d epends on 
the mass accretion rate and black hole mass (e.g., iGoodmanl 
120031: iKing et al.H2008l) . In what follows, for illustration pur- 
poses, we adopt fiducial values of a = 0. 1, h = 10~ 2 , f m = 10~ 2 , 
R out = 10 4 # g , M. = 1O 8 M , and a = 0.998, where R g = GM./c 2 
is the gravitational radius. 

3.1. Transporting Angular Momentum 

In planar accretion disks, matter gradually falls inward 
the central black hole, whereas nearly all of the angu- 
lar momentum is transported outward due to viscosities 
(lFranketal.1119921). Thi s remains valid for warped accretion 
disks (IScheuer & Feilerlll996b . The difference is that gravit- 
omagnetic interactions of the hole and the disk induce extra 
transfers of angular momentum between them. This transfer 
rate depends on the Lense-Thirring precessing rate and warp 
propagation velocity. There exists a transition radius where 
the Lense-Thirring precessing timescale is comparable with 
the warp propagation timescale, which is usually denoted by 
a "warp radius". 

Within the warp radius, the Lense-Thirring precessions are 
dominant. Viscosities damp out the differential precessions 
rapidly, giving rise to an inner flat disk. It is easy to verify: if 
initially the angle of the disk and the hole is smaller than ir/2, 
the inner flat disk is aligned with the hole. By contrast, if 
the angle is larger than ir/2, the inner flat disk is anti-aligned. 
In the region around the warp radius, the incoming matter is 
inclined toward the inner flat disk, and the associated angular 
momentum is transferred to the hole through gravitomagnetic 
interaction. Further beyond the warp radius, warps propagate 
outward efficiently and the misaligned angular momentum is 
transported outward as well. 

Since the gravitomagnetic torque integrated over the whole 
disk in Equation © is always orthogon al to Jh, it can be ex- 
pressed in the form of (King et al. 2 005|) 

T = -2tt y T LT RdR = -K L J h x J d - K 2 J h x (J h x J d ), (7) 



where a hat symbol on top of a vector stands for the corre- 
sponding direction vector, K\ and K 2 are time-dependent co- 
efficients, and 



- 2 "i 



LRdR, 



(8) 



is the total angular momentum of the disk. Here, the first term 
corresponds to the precessions and the second term corre- 
sponds to changes in the angle between the hole and the disk. 
Below we will numericall y demonstrate that K\ and K 2 are 
generally positive ( see also fecheuer & Feilerll9 96: Kin g et all 
2005; Martin et al. 2007). 



3.2. Magnitude Estimates 

In in previous section, we demonstrate that the charac- 
teristic extension of warps is determined by the timescales 
for warp propagation given by Equation © and for the lo- 
cal Lense-Thirring prec essing tur ~ |L|/|T LT | = c 2 R 3 /2GJh 
(IScheuer & Feilerlll996l) . By equating these two timescales, 
one obtains the warp radius R w = 2GJ ri /v2C 2 . Using the for- 
mulation of V2 and the fiducial values of parameters, R w is of 
the order of 



R« 



V h 2 J 



2/3 



10- 



-4/3 



(9) 



Considering that the Lense-Thirring torque in Equation (2} 
falls off rapidly with the radius in proportion to R~ 3 , the major 
contribution it has on the hole comes from the region of the 
disk around R w . As a result, the torque that the disk exerts on 
the hole approximates 

4irGJ h L(R w ) 



Rw 



V2 



M^GM.R^ 



(10) 



where the mass accretion rate M ~ ?>kv\ E is used for simplic- 
ity. This results in a timescale fo r the alignment bet ween the 
hole and the whole disk (see also lPerego et al.ll2009|) 

J h V\M. [Ry 

h\ = — ~a ~y\ 

Tlt vi M V R 



io-W gr (^ 




(ii) 



where t gY = M./M is the growth timescale for a black hole 
accreting at a rate of M. Again, with the help of the approxi- 
mation M ~ 37r^iE, one can estimate the mass accretion rate 
as 

af m M.h 2 



Mr 



R g (R out /R g ) 3 / 2 
h 



10" 8 M. yr" 



\0.lJ 



10" 



fm \ ( _^o 



10- 2 J V 10% 



-3/2 



(12) 



and the growth timescale t gY ~ 10 8 yr. We note that this esti- 
mated accretion rate is generally equal to the Eddington limit. 

If there is no gas supply to the accretion disk, the disk will 
be consumed on a timescale of an order of 



.Md 
M 



fmM. 

M 



10 6 yr 



(id' 



10" 



V 10% 



3/2 
(1?) 



which is independent of f m and M m . As can be seen from 
Equation (fTTT) . once the disk is gradually depleted (without 
gas replenishment), the resulting decreases of surface den- 
sity and the mass accretion rate will significantly prolong 
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the alignment processes. The timescale for warp propagation 
throughout the disk in Equation © is 



t w =^ = 2ah- 2 n~ l = 3x 10 4 yr 



a \ ( h 



(£) 



( M. \ /_R, 



io- 2 y V 1Q8M / V 1Q4 ^g 

The total angular momentum of the disk is 



3/2 



• (14) 



Jd = |Jd| ~ ■jfmM.y/GM.Rwt. (15) 



This gives a ratio of Jd over /h- 



3 /m / ^out _ 3 i ( fm \ ( _^ 



J h 4 a\l R g C V 10-2 ) \ W 4 R 



1/2 



(16) 



3.3. Alignments and Anti-alignments 

In this section, we qualitatively study what controls the fi- 
nal configuration of black hole-disk systems, i.e., alignments 
or ant i-alignments, inspired by the arguments of iKing et ail 
(|2005|) . As shown in the previous section, the (anti-)alignment 
timescale is generally orders of magnitude less than the 
growth timescale for black holes. Therefore, the accretion 
contribution to the changes of Jh in Equation © can be safely 
neglected. By denoting mc as the angle between Jh and Jd 
and using Equation ©, we have 



dt 



( Jh • Jd) = Jh • + K 2 J d sin 2 iT 
at 



(17) 



In the meantime, above the time derivative can be alterna- 
tively written 

d dJ d d 

-r ( Jh • Jd) = Jh-r- cos mc + J h J d — (cos inc ) , (18) 

dt dt dt 

where we apply the fact that the gravitomagnetic interaction 
only modifies the direction of Jd but keeps its modulus un- 
changed (i.e., dJh/dt = 0). A simple algebra manipulation 
yields 

d 2 dJd d/d 

AJd -7- (cos inc ) = K 2 J d sin inc + J h • — - - J h — cos inc . 
dt dt dt 

(19) 

Obviously, besides the gravitomagnetic interaction, the evolu- 
tion of the inclination angle mc depends on how the angular 
momentum of the disk Jd changes. 

According to the different types of disk feeding in the outer 
edges as illustrated in Figure Q] we distinguish the two follow- 
ing cases. 

1) Fixed flows. If the disk is steadily fed and maintains a 
preferred angular momentum distribution, one shall expect Jd, 
and therefore J d , to remain somewhat unchanged with time. 
As a result, Equation (TT91) is simplified into 



(cos6> inc ) = ^ sin 2 6> inc > 0. 
dt 7h 



(20) 



Thi s means that, as found by previous in t ensive stud- 
ies (iScheuer & Feileri[l99l iMartin et al.l 120071: fPerego et al.l 
120091) . a continuously fed disk always drives the hole to align 
with it regardless of the initial inclination. 

2) Free flows. For disks without replenishment, the total 
angular momentum of the system (Jh+ Jd) is conserved. With 



the help of Equa tion ([7]), it is eas y to derive the following 
equation (see also lKing et al.ll2005h 



d ^ 2 sin h 

— (cos6> inc )= 

dt Jh 



l + ^COS0inc 
•Al 



(21) 



Note that here the coefficient K 2 is time-dependent. Because 
K 2 is positive, once initially cos#i nc > 0, i.e., 6 mc < ir/2, the 
right hand side in the above equation is never smaller than 
zero. Therefore, 9 mc will decrease continuously and the sys- 
tem ends up aligned. For #i nc > 7r/2, there exist two subcases: 

• if Jd/Jh > - cosine, the inclination mc decreases; 

• while on the contrary if J d /A < -cos#i nc , the inclina- 
tion mc increases toward tt. 

In real situations, J d generally decreases with time due to the 
internal viscous dissipation, and subsequently, the ratio Ja/A 
can transit from the former subcase to the latter. Therefore, 
the above conditions are indeed instantaneous for the behavior 
of the inclination. We next dig into exploring the detailed 
alignment/anti- alignment processes with help from numerical 
calculations. 

4. NUMERICAL SCHEME 

Numerically solving Equations (Q]) and © gives the self- 
consistent evolution of the black hole accretion disk system. 
For this purpose, w e implement a differencing scheme follow- 
ing jPringiJ ([1992), but with some different treatment tricks. 
A detailed description of our differencing scheme and robust 
tests of the code are given in Appendices A and B, respec- 
tively. 

In order to produce the free-torque conditions in the in- 
ner boundary, we enforce L,(Rm) = such that the mass 
and angular momentum that reaches the inner boundary 
are removed, resulting in accretion onto the hole (see also 
iBregman & Alexander! I20II iNixon & Kingll2012l) . Accord- 
ingly, the initial surface density is setup by ad ding an ex- 
tra fa ctor to allow for the torque-free condition dFrank et al.l 
fl992h . namely, 



X(R) = X R- p 




^R^Ra 



(22) 



where R[ n is the inner edge of the disk. In consideration of the 
inner portion of the disk aligned/anti- aligned with the hole, 
we set di/dR = 0dXR m . 

Given an initial inclination angle of #i nc ,o, we generate the 
initial distribution of the direction of disk rings as 

i(R) = (sin 

^inc,0? 0, COS $inc,o)? f ° r ^in R ^out- (23) 

For this configuration, we presume that initially the black hole 
spin is oriented toward the z-axis. 

The outer boundary conditions at R out are treated differently 
for fixed flows and free flows. 

• Fixed flows. Gas supply at the outer edge inherits a 
specified direction of angular momentum, that is to 
say £(R 0U t) keeps unchanged throughout. The surface 
density at the outer boundary is fixed to T>(R out ) = 

Eo^out ( l ~ V^in/^out) • This yields a mass supply 
through the outer boundary at a rate given by Equa- 
tion dl2t . 
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Fig. 2. — Shape of accretion disks with initial inclinations of #i nc ,o : 
viscous timescale at R on t (see Equation ©). 



30° (top) and 150° (bottom) for free flows at various times. Here t v = 2 x 10 6 yr is the 



• Freeflows. Angular momentum freely diffuses outward 
through R out . To ensure a correct treatment of the dif- 
fusion at R 0VLt , we choose the outer grid boundary well 
beyond the outer disk edge (e.g., ~ 1007? out ). 

In our calculations, without specified otherwise, we adopt 
the fiducial values a = 0.1, p=l/2,h= 1(T 2 , R out = 10 4 # g , 
Ri n = 6R g , and the initial black hole mass M. = 1O 8 M and 
spin a = 0.998. In principle, R{ n is determined by the marginal 
stable orbit of the central black hole. However, since the 
inner portion of the disk is always aligned/anti- aligned with 
the hole, the location of R[ n does not affect the gravitomag- 
netic torque and hence the alignment rate. Moreover, a rel- 
atively larger R{ n helps to improve the time consumption of 
the computations in view of the timestep size limited by vis- 
cous timescales (see Equation ©). To assess the processes 
of alignments and anti- alignments, we adjust the two free pa- 
rameters: the mass fraction of the disk to the black hole f m , 
which determines the surface density of the disk in Equation 
(l22l) and hence the mass accretion; and the initial inclination 
angle between the disk and the hole 0i nc ,o- We setup the grid 
with 202 points and evolve the equation set for a time of \.5t v , 
where t u = 2x 10 6 yr is the viscous time at R out (see Appendix 
A for details). 

5. RESULTS 

5.1. Evolution of Warped Accretion Disks 

We define the inclination of disks at radius R with respect 
to black holes as 

_i Jh-L(fl) 



In Figure [2 we illustrate how the shape of warped accretion 
disks evolves with times for free flows. The initial inclina- 
tion angles are #i nc ,o = 30° in the upper panel and ft nCj o = 150° 
in the bottom panel, respectively. We can find that an inner 
flat disk rapidly forms for all cases and is aligned to the black 
hole for #i nC5 o = 30° and anti-aligned for #i nC5 o = 150°. The 
size of the inner flat disk continuously grows until roughly 
l0 2 Rg, corresponding to the warp radius R w defined in Equa- 
tion ©. Warp propagation reaches the outer edge of the disk 
within a timescale of ~ 0.01^, consistent with the estimate 
from Equation (fT4l) . The disk for 6 mc $ = 30° gradually ap- 
proaches full alignments; however, the alignment rate is prone 
to severe attenuation with time. The disk inclination is almost 
unchanged from 1 .0 to l.5t y . At the end of the numerical cal- 
culation 1.5t u , the disk maintains an inclination of mc ~ 2° 
at the outer edge. Henceforth, it hardly aligns with the black 
hole fully. This is because of a severe depletion of the surface 
density due to mass accretion (see Equation ([T3])). In Fig- 
ure it is verified that the mass accretion onto the black hole 
decreases rapidly with time by more than one order of mag- 
nitude as a result of disk depletion. This implies that if there 
is no gas supply, the luminosity of the accretion disk will fade 
out less than 10 6 yr. We further calculate the components of 
the gravitomagnetic torque that the whole disk exerts on the 
hole K\ and K2 as defined in Equation ©. As shown in Figure 
IU it is confirmed that irrespective of the initial inclinations, 



both K\ and K2 are positi ve. Moreover, K\ and K2 are in the 
same order of magnitude 
l200llMartIn^tl^ 



mc (R) = cos 



|Jh||L(tf)| 



(24) 



The system for 6^i nC; o = 150° evolves toward anti-alignment 
in a similar fashion, as expected since initially Ja/Jh = 0-78 < 
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1 2 3 

Time (10 8 yr) 

Fig. 3. — Mass accretion rate onto black holes (M m ) with time for free 
flows. 

cos6>i nC5 o = 0.87. However, the inclination 6 mc at R out has 
changed only from 150° to ~ 160° at the end of the cal- 
culation, which is very inefficient compared to the case of 
#inc,o = 30°. The corresponding reason is that, as shown in 
Equation (BTT) . the change of mc from the gravitomagnetic 
torque is proportional to ex (1 +Jh/Jd cos#i nc ). Its value for 
#inc,o = 150° is smaller by an order of magnitude than that for 

#inc,0 = 30°. 

The cases for fixed flows are illustrated in Figure [5] By 
analogy to free flows, warps spread to the outer boundary 
within 0.01^; however, the subsequent evolution is some- 
what different. The system with 6{ nC: o = 30° quickly achieves 
full alignment after a time less than l.0t u . For #i nc? o = 150°, 
the disk shape shows complicated behaviors. The inner anti- 
aligned angular momentum is transported outward by warp 
propagation and confronts the angular momentum carried by 
the feeding matter. As a result, an abrupt discontinuity of the 
inclination appears and the disk breaks into two parts: an in- 
ner nearly ant i-aligned portion and an outer misaligned por- 
tion (see also lLodato & Pringlel 120061: iNixon & King 120121 : 
iNixon et all 120 121) . The location of the discontinuity trans- 
fers inward with time and the inner portion will be eventually 
swallowed by the hole due to mass accretion. In the mean- 
time, the black hole is driven to align with the outer portion 
progressively. Around the time 0.34-0.38^, after the in- 
ner anti-aligned portion of the disk is consumed, the newly 
formed inner disk turns to be aligned with the black hole. At 
the end of the calculation 1.5t u , there is complete alignment 
between the disk and the black hole. 

The mass accretion rates onto the black holes for fixed flows 
are plotted in Figure|6] From Equation (|2]>, the role of warping 
increases the inflow velocity and therefore enhances the mass 
accretion rate. As expected, for #i nC5 o = 30°, the accretion rate 
onto the hole (Mm) has a mildly decreasing trend by a factor of 
two at the beginning in response to the approaching alignment 
of the system. Then it maintains a steady value of ~ lM yr _1 , 
equal to the gas supply rate at the outer boundary. Quite dif- 
ferently, for = 150° there exist two peaks of the mass 
accretion rate at times ~ 0.2^ and ~ 0.5^ due to the large 
magnitude of warping developed at these moments (see the 
bottom panel of Figure [5]) . While the inner anti-aligned disk 
portion is being consumed, the accretion rate becomes signifi- 
cantly attenuated by orders of magnitude and reaches the min- 
imum (Mi n — >> 0) around 0.34-0.38^. Nevertheless, the mass 
accretion onto the hole eventually approaches a steady rate, 




Time (10 6 yr) 

Fig. 4. — Components of the gravitomagnetic torque in Equation J7} for 
free flows. K\ corresponds to the precessions and K2 corresponds to aligning 
the system. Note that both K\ and K2 are positive. 

the same as for 6{ nC: o = 30°. It is worth pointing out that the 
black hole undergoes retrograde accretion for a time ~ 0.35^, 
after which it transits to prograde accretion. 

5.2. Alignment and Anti- alignment Time scales 

Figure |7] plots the evolution of the inclination angles be- 
tween the black hole and the disk for two classes of accretion 
flows. Here the inclination angle is defined by 



= cos 



Jd • Jh 

iJdljjhT 



(25) 



where Jd is the total angular momentum of the disk. The 
upper left panel shows mc /Q mc $ against time with 6{ nC} o = 
(10°, 30°, 40°, 50°, 60°) from bottom to top for free flows. 
The timescales that the inclination angle decreases to 1 je of 
its initial value generally lie at several 10 5 yr, consistent with 
the magnitude estimate in Equation (TTTb . However, the decay- 
ing rate of the inclination angle rapidly slows down with time 
because of the depletion of the disk (almost all of the disk 
angular momentum is transported to the large radius). Af- 
ter ~ 10 6 yr, the inclinations for all the cases are prone to be 
frozen, indicating that the disk shall hardly approach complete 
alignment. This further confirms the results in the previous 
section. The upper right panel plots (n - 0{ nc ) /(tt- 0{ nc 0) with 
0inc,o = (120°, 130°, 140°, 150°, 170°) from top to bottom for 
free flows. Overall, although quite inefficient, the inclination 
angles have a tendency to approach 8 mc —> 180° except for 
#inc,o = 120°, implying a trend of anti- alignment between the 
hole and the disk. The system for 6^i nC; o = 120° tends toward 
alignment instead of anti-alignment because Jd/A> -cos#i nc 
is always retained throughout the calculation (see Section 
3.3). Interestingly, there is an increasing trend at the begin- 
ning for #i nC5 o = 130° and 140°, because Jd/A> -cos#i nc ini- 
tially, which, however, transits to Ja/A < -cos^ nc later on. 
Again, the alignment rates are significantly slower compared 
with the counterparts in the upper left panel. 

The bottom panels of Figure |7] show the results for fixed 
flows. The timescales where mc decreases to l/e of its initial 
value are also generally several 10 5 yr for Q mCi o < 90° (bottom 
left panel) and almost insensitive to #i nc ,o- Within a time of 
1.0^, the entire disks have been in complete alignment with 
the holes. In the bottom right panel, mc decreases with time 
even though Q mc $ > 90° as was expected through the previous 
theoretical analysis. There are inflection points in the evolu- 
tion curves of mc that correspond to when the disks break into 
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Time (10 6 yr) 

Fig. 6. — Mass accretion rate onto black holes (M- m ) with time for fixed 
flows. The tiny ripples in the 0- mc ^ = 150° line is due to the numerical effect. 

two parts as illustrated in Figure [5] Generally, the timescales 
of the decaying rate for 6{ nC} o > 90° (bottom right panel) are 
relatively longer compared to these for < 90° (bottom 
left panel) due to the same reasons for free flows. However, 
because the disks are replenished and therefore the gravito- 
magnetic torques from the disks are not attenuated, the de- 
creasing rates of mc are clearly more rapid than their coun- 
terparts (top right panel) for free flows. 

In summary, for fixed flows in which disks are continually 
fed, alignments always occur regardless of the initial inclina- 
tion. The alignment timescale is of an order of 2 x 10 6 yr. For 
free flows in which there is no gas supply, both alignments 
and anti-alignments could occur. However, as a result of the 
notable depletion of the disk's surface density at a timescale 



of 10 6 yr, the alignment/anti- alignment rate is significantly re- 
duced, leading the disks to be not completely aligned/anti- 
aligned. 

5.3. Conditions for Anti-alignments 

Previous sections show that anti- alignment between black 
holes and disks is possible only when the initial inclination 
angle is 6 mc $ > 90° for free flows. To assess how the final 
configuration of the system depends on the free parameters, in 
this section we evolve the equation set given different param- 
eters f m and #i nC5 o only for free flows. We note that the param- 
eter f m determines the initial total angular momentum of the 
disks as in Equation (IT6l) . Figure [8]plots the evolution of incli- 
nation angle 9 mc between Jh and Jd with f m /\0~ 2 =0.1, 1, 1.5, 
2, 3, 4, and 6, corresponding to the initial ratios Jd/A = 0.078, 
0.78, 1.17, 1.56, 2.34, and 3.12, respectively. The initial in- 
clination angles are #i nC5 o = 160°, 140°, and 120° from top 
to the bottom panel. Since gas diffuses outward through R out 
for free flows, there is an outflow of angular momentum at 
R out accordingly. Here Jd includes these aspects of angular 
momentum. For comparison purposes, we also calculate the 
inclination angles in terms of the angular momentum of the 
disk just within R ouU plotted by dotted lines in Figure [U As 
can be seen, the evolution trends of mc depend both on the 
ratios Ja/Jh and the initial inclinations. Specifically, there are 
three prominent features as follows. 

1) Consistent with the prediction in Section 3.3, the in- 
clination angles 6 mc generally go up when the initial ratios 
Jd/Jh < -cos#i nC5 o, t> ut otherwise, drop off. We note that in 
the middle panel of Figure [8] the inclination of the system 
with an initial ratio Ja/A = 0.78, larger than -cos 6^i nC; o = 0.76, 
descends mildly at the very beginning and then increases to- 
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Fig. 7. — Evolution of the inclination between the holes and the disks for two classes of boundary conditions: (upper) free flows and (bottom) fixed flows. Two 
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Fig. 8. — Evolution of the inclination angle bet ween the angular momenta of the black holes and the disks for free flows with different initial Jd/Jh, which 
is determined by the parameter f m as in Equation ( fl6l The initial inclination angles are 6- m c,o = 160° (top), 140° (middle), and 120° (bottom). Dotted lines 
correspond to the angles in terms of the angular momentum of the disk just within R ou t (excluding the parts that outflow through Rout). Gray horizontal dashed 
lines correspond to mc = 90° . 



ward 180°. This is because the ratio Jd/Jh rapidly decreases 
to Jd/Jh < -cos# inc . 

2) The inclination for ft nC5 o = 160° begins to enter 6 mc < 90° 
when the initial ratio Jd/Jh is larger than ~ 3. This means that 
the system will approach alignment, and most importantly, the 
black holes undergo a transition from retrograde to prograde 
accretion. Similarly, the critical ratios for 6 mc $ = 140° and 
120° are ~ 2.0 and 1.2, respectively. 

3) Without replenishment, the disks are greatly depleted on 
a timescale of 10 6 yr. As a result, the gravitomagnetic torques 
are reduced with time and the inclination angles are almost 
frozen to values not equal to 0° or 1 80° at the end of calcu- 



lations. In other words, the systems hardly achieve full align- 
ment or anti-alignment. 

On the basis of the angular momentum conservation of the 
system (free flows in the present study), iKing et ail (120051) 
proposed that anti- alignments will occur provided that the ini- 
tial angle between the angular momenta of the hole and the 
disk satisfy Jd/Jh ^ -2cos#i nC5 o. Applying this formula to the 
present cases, for Q- mc $ = 160°, anti-alignment occurs when 
the ratio obeys Jd/Jh < 1.88 initially; for Q mCi o = 140° and 
120°, the corresponding critical ratios are 1.53 and 1.0, re- 
spectively. Our results in Figure [8] are in considerable dis- 
agreement with these values. We ascribe such discrepancy 
to the neglect of the disk's depletion in the previous work. 
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The ratio between the timescales for (anti-)alignments and 
the disk's depletion is of an order of unity for fiducial pa- 
rameters from Figure |7] During accretion onto the black hole, 
the disk's mass is gradually consumed and the majority of the 
disk's angular mo mentum is transported to the large radius 
dFrank et al.lll992|) . Note that the gravitomagnetic torque is in 
proportion to 7T 3 , therefore the (anti-)alignment rate will be 
greatly reduced. This effect leads to some systems with ratios 
of Jd/A > -2cos#i nC5 o still maintaining 9 mc > 90° at the end 
of the calculations in Figure [U 

6. DISCUSSIONS 
6. 1 . Implications for Black Hole Spin 

In addition to fueling black holes, accretion also adds an- 
gular momentum carried by matter to black holes and ac- 
cordingly modifies their spin. How spin changes depends on 
the orbits at the inner edge of the accretion disks. Specif- 
ically, if the disk is aligned with the black hole, prograde 
accretion will spin up the hole; otherwise, retrograde accre- 
tion spin down the hole. The black hole will become max- 
imally rotating (a —> 1) when its mass is doubled provided 
prograde or alternative retrograde acc retion is prefe rably re- 
tained for a sufficient amount of time (Thorne 1974]). In view 
of the episodic and random activi ties of black hole accretion 
(IWang et al.l 120061: iLi et al.ll20loh . it is crucial to model the 
spin evolution whether the accretion in each episode proceeds 
progradely or retrogradely. 

Our results obtained above imply that the configurations 
of the accretion disks depend on feeding them at the outer 
boundaries. If the disks are continually fed for a long enough 
lifetime (say, e.g., 10 7 yr), alignments always occur. In this 
case, prograde accretion is the dominated fashion. We ex- 
pect that after a series of activities, black holes will be close 
to maximal rotation. On the other hand, if there is no gas 
feeding, both alignments and anti- alignments are possible, de- 
pending on the initial inclinations and the ratios Jd/A (see 
Section 5.3). One has to follow the probabilities of align- 
ments and anti-alignments over episodes using the method 
presented here to determine spin evolution. An incorpora- 
tion of the present calculations into semi-analytical models 
of galaxy formation and ev olution will shed light on the exist- 
ing studies in the field (e.g..lLagos et a l. 2009; Fanidaki s et al.l 
l2011l:lBaraussell2012l:rvblonteri et al.ll20ll iDotti et al.ll2013b . 
and provide more realistic information on SMBH spin. 

Note that black hole accretion only lasts a lifetime of ~ 
10 6 yr as a result of disk depletion without feeding. 3 Indeed, 
the lifetime of feeding disks directly determines the observ- 
able episodic lifetimes of AGNs. In this sense, observational 
measurements of the episodic lifetimes of AGNs can place 
useful constraints on the SMBH spin evolution. Unfortu- 
nately, thus far the episodic lifetime of AGNs remains elu- 
si ye though som e endeavor has been put forth (see a review 
of iMartiriil 12004)) . The most compelling method — proximity 
effects estimated the episodic lifetimes in a broad range of 
~ 1 6 - 10 8 yr, varying from study to study dKirkman & Tvtlerl 
120081: iGoncalves et al.ll2008tlFurlanetto & Lidzll2011h . Future 
observations from large sample quasar surveys with a vari- 
ety of techniques would offer promising measurements on the 

3 The lifetime can be generally estimated as follows: for black holes ac- 
creting at the Eddington limit MEdd ~ IM^Mqjy' 1 , where Mg = M / 1O 8 M0 , 
the accretion disks with masses of ~ (H/R)M. = 10 6 M 8 M© will be con- 
sumed at ~ 10 6 yr, where H/R^10~ 2 . 



episodic lifetimes and help to understand the cosmological 
evolution of SMBH spin. 

On the other hand, we turn over the issue so as to con- 
strain the episodic lifetimes of AGNs from the otherwise ob- 
tained SMBH spin evolution. Recently, quantifying the radia- 
tive efficiency of mass accretion through SMBH demography 
strongly suggests tha t SMBHs are spinning down with cos- 
mic time since z ~ 2 dWans et al. 2009; Li et al. 20 EE l20H 
see also the discussion of IZhang et al. 201 2). Combined with 
the present results, this indicates that the episodic lifetime of 
AGNs must not exceed several 10 6 yr and probably have cos- 
mological evolution as well. More detailed modelings on the 
connections between SMBH spin evolution and the episodic 
lifetime of AGNs are highly deserved in future works. 

It is worth stressing that we presume black holes to be 
accreting at the Eddington limit in the calculations. The 
timescale for (anti-)alignments between the disks and the 
holes is in proportion to the reciprocal of the mass accretion 
rate (see Equation (fTTT)). Therefore, if the presumed mass ac- 
cretion rate is relaxed, the above proposed limit on the AGN 
lifetime changes in the same manner accordingly. 

6.2. Uncertainties 

In our calculations, the inner edge of the disk is fixed at 
R in = 6R g . Considering that the inner portion of the disk within 
R w always remains aligned with the hole, this does not affect 
the gravitomagnetic torque on the hole and therefore the align- 
ment rate. On the other hand, during the alignment course, 
the accreted mass and angular momentum are negligible with 
those of the black hole. Hence, the results presented in this 
paper are insensitive to the location of the inner boundary. 
For the sake of simplicity, the disk's outer edge (limited by 
the self-gravitating) is fixed at R out = \0 4 R g . In reality, the 
self-gravitating radius depends on the accretion rate, but gen- 
erally lies at a range of ~ 10 3 - 10 4 ^ g for stan dard accretion 
disks (e.g.. iGoodmanl 120031: lUng et al.ll2008h . The location 

of R out just determines the ratio of Jd/Aoc R l J^ as in Equa- 
tion (TT6K but does not affect the gravitomagnetic torque and 
(anti-)alignment timescale as in Equations © and (fill) . The 
overall results remain unchanged. 

The viscosities for accretion disks are treate d using the first- 
order approximations given by lO gilviel (1 1 999b . The nonlinear 
fluid effects that appear for large amplitude wa rp would mod- 
ify the expressions of the effective viscosities (O gilvielfl999b 
lLodato & Pried l20lo|) . We expect this to become significant 
when the disks break as shown in the bottom panels of Fig- 
ure[5] In this case, however, it is unclear whether the evolution 
equation set adopted here remains adequate. Dedicated inves- 
tigations, including nonlinear fluid dynamics, are beyond the 
scope of the present study. We defer this to future sophisti- 
cated numerical simulations. 

Finally, we consider the accretion disks under the Newto- 
nian timespace and treat the gravitomagnetic interaction in the 
weak-field limit (the Lense-Thirring precession). The general 
relativistic calculations are of somewhat less importance be- 
cause the warp radius R w is far beyond the gravitational radius 
R g (see Equation ©). 

7. CONCLUSIONS 

We numerically solve the evolution equation set of spinning 
black holes and their warped accretion disks. To mimic the 
realistic accretion processes, we consider two classes of ac- 
cretion disks in terms of the outer boundary conditions: fixed 
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flows for which the disk is continually fed with a preferred 
angular momentum through its outer edge; and free flows for 
which there is no gas supply and the disk diffuses freely at its 
outer edge (see Figure [TJ. Our main results are as follows. 

1 . For disks continually fed, alignments between the holes 
and disks always occur at a time of several 10 5 yr re- 
gardless of the initial inclinations. If the initial inclina- 
tion angles are #i nC5 o > tt/2, black hole accretion transits 
from retrograde fashion to prograde fashion after a time 
of several alignment timescales. The mass growths dur- 
ing these two phases are comparable. 

2. For disks without a gas supply, both alignments and 
anti-alignments are possible, depending on the initial 
inclinations and the ratios Ja/A- It is worthwhile 
to point out that the lifetime of black hole accretion at 
the Eddington limit in such cases last ~ 10 6 yr due to 
the depletion of the accretion disks. This has important 
influences on the (anti-)alignment processes. 



3. Black hole spin is closely linked to the manners of mass 
accretion, namely, progradely or retrogradely. Care- 
ful inclusion of disk feeding at the outer boundaries 
is highly important to modeling spin evolution in the 
context of episodic and random activities of SMBHs. 
The episodic lifetime of AGNs, which can be measured 
with a variety of techniques, is directly determined by 
the lifetime of disk feeding. We therefore propose that 
there would be a close connection between the black 
hole spin and the episodic lifetime of AGNs, which de- 
serves a comprehensive investigation in future works. 
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APPENDIX 
A. DIFFERENCING THE EQUATION SET 

We construct the differencing scheme of Equations (OQ) and © generally following iPringlel (1 1 9921) . but with some modifications. 
For the sake of completeness, we summarize our implementation here. 

The spatial domain is discretized into (1+2) logarithmic grid points: Rf = R{ n e lAx , for i = 0, ...,/+ 1, with Ax the equal spacing 
of the logarithmic grid. The points Ro and R\+\ represent the inner and outer grid boundaries, respectively. If we define the 
advective velocity 



dt 



dR 



and use Equation (|2]), Equation (Q]) can be written as 

dL 

~di 



Id Id 

--^(Vadv*L) + - — 



d 1 di 

3^-(^ + -^L- 



ld_ 
RdR 



V3RL x 



di 
dR 



2G J h x L 



(Al) 



(A2) 



Notationally, we use a su perscript n to denote the timest ep and a subscript / to denote the spatial grid point. The differencing 
scheme of Equation (IA2b is then built as follows (see also Bre gman & Alex ander 20T2I) 



n+l . 



(VW&.JfcrtLk -OWW , 3 { - C 1/2 - [W-^du] C 1/2 } 



At 



RjAx 



Rf(Ax) 2 



2R 2 (Ax) 2 



Rf(Ax) 2 



c 2 R? ' 



(A3) 



where At is the timestep size, and the advective term is treated using upstream differencing: for V ac i v >0,k = i—l; whereas for 
V a dv < 0, k = i. The advective velocity is calculated by 



(Vadv)? = ^ 



Similarly, the differencing scheme of Equation © is built as 



W+i 



W-l 



2Ax 



(A4) 



Jn+l _ jn 
_h 

A* 



4ttG 



L7 



= (Minjin) n -JExV^Al 



(A5) 



where the mass accretion rate onto the hole is calculated as (M- m ) n = -27t(7?VrE)" and the specific angular momentum at the 
inner edge is calculated as jf n = (R 2 ^})"^. As described in Section 4, at the inner boundary, we take Lq = 0, and £q = t x \ at the 
outer boundary, for free flows, we take L" +1 = (z/iL)"/(z/i)"; for fixed flows we take L" +1 = ^ +1 fi/+iEo/?7+i(l _ y/Rm/Ri+i) and fix 
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Fig. 9. — To verify the validity of our numerical code, three tests are performed (see the text for details). Left: the surface density of a flat disk (#i nc ,o = or tt) 
at the beginning and end of the calculations. Middle: inclination profile of the disk in a steady state with the outer inclination angle fixed at mc = 10°. Right: 
Evolution of inclination angle between the hole and the disk with an initial angle of 0- mc fi = 10°. 

t M = (sin#i nC5 o,0,cos#i nC5 o). This enforces a coherent mass supply at a rate given roughly by Equation (fT2l) . The timestep size is 
adjusted for every step according to 



where 



At = - min {Ah , At 2 , At 3 , A* C fl , At LT , At h } , 



. (RjAxf . 
At; = min — - — - — lor / = 1,2,3, 

J !<«</ OA- 



(A6) 



(A7) 



A* 



cfl = mm 



RiAx 

1<7<7 (Vadv)/ ' 



At] 



lt = min 



i<«</ T LT ; ' 



and Afh = A 



4ttG t y-L, a 

Z=l 



(A8) 



B . VALIDITY OF THE NUMERICAL CODE 

To verify the validity of our numerical code, we perform three tests. (1) We run the code upon a flat disk (0in C ,o = or 7r) for a 
time of 1.5^. The left panel of Figure [9] plots the surface density of the disk at the beginning and end of the calculations. There 
are tiny deviations less than 2% due to the dissipation of the numerical scheme. (2) For small warping, the re exist analytical 
solutions if disks are continually fed at the outer edges, analogous to fixed flows defined in the present paper (IScheuer & Feilerl 
1996; M artin et aLll2007l: [Chen e t al. 2009|). We only evolve Equation (OQ) with a sufficient amount of time using the differencing 
scheme described above and obtain a steady shape of the disk. In the middle panel of Figure we compare th e numerically 
obtained inclination profile of the disk with the analytic solution given by Equation (24) in iMartin et al.l (|2007). A complete 
match can be found. (3) Then initializing the disk with the above shape, we go on to evolve Equations (Q]) and © simultaneously 
to obtain the time-dependent inclination angle between the hole and the disk. The right pan el of Figure[9]shows that our numerical 
calculations are in exact agreement with the analytical solutions given by Equation (52) in lMartin et all(l2007l) . These three tests 
indicate a validity of our numerical code. 
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